Search for: All records

In this paper, we derive correlations between core-collapse supernova observables and progenitor core structures that emerge from our suite of 20 state-of-the-art 3D core-collapse supernova simulations carried to late times. This is the largest such collection of 3D supernova models ever generated and allows one to witness and derive testable patterns that might otherwise be obscured when studying one or a few models in isolation. From this panoramic perspective, we have discovered correlations between explosion energy, neutron star gravitational birth masses,⁵⁶Ni andα-rich freezeout yields, and pulsar kicks and theoretically important correlations with the compactness parameter of progenitor structure. We find a correlation between explosion energy and progenitor mantle binding energy, suggesting that such explosions are self-regulating. We also find a testable correlation between explosion energy and measures of explosion asymmetry, such as the ejecta energy and mass dipoles. While the correlations between two observables are roughly independent of the progenitor zero-age main-sequence (ZAMS) mass, the many correlations we derive with compactness cannot unambiguously be tied to a particular progenitor ZAMS mass. This relationship depends on the compactness/ZAMS mass mapping associated with the massive star progenitor models employed. Therefore, our derived correlations between compactness and observables may be more robust than with ZAMS mass but can nevertheless be used in the future once massive star modeling has converged.

We study in detail the ejecta conditions and theoretical nucleosynthetic results for 18 three-dimensional core-collapse supernova (CCSN) simulations done by Fornax. Most of the simulations are carried out to at least 3 s after bounce, which allows us to follow their longer-term behaviors. We find that multidimensional effects introduce many complexities into the ejecta conditions. We see a stochastic electron fraction evolution, complex peak temperature distributions and histories, and long-tail distributions of the time spent within nucleosynthetic temperature ranges. These all lead to substantial variation in CCSN nucleosynthetic yields and differences from 1D results. We discuss the production of lighterα-nuclei, radioactive isotopes, heavier elements, and a few isotopes of special interest. Comparing pre-CCSN and CCSN contributions, we find that a significant fraction of elements between roughly Si and Ge are generally produced in CCSNe. We find that⁴⁴Ti exhibits an extended production timescale as compared to⁵⁶Ni, which may explain its different distribution and higher than previously predicted abundances in supernova remnants such as Cas A and SN1987A. We also discuss the morphology of the ejected elements. This study highlights the high-level diversity of ejecta conditions and nucleosynthetic results in 3D CCSN simulations and emphasizes the need for additional long-term (∼10 s) 3D simulations to properly address such complexities.

Using 20 long-term 3D core-collapse supernova simulations, we find that lower compactness progenitors that explode quasi-spherically due to the short delay to explosion experience smaller neutron star recoil kicks in the ∼100−200 km s⁻¹range, while higher compactness progenitors that explode later and more aspherically leave neutron stars with kicks in the ∼300−1000 km s⁻¹range. In addition, we find that these two classes are correlated with the gravitational mass of the neutron star. This correlation suggests that the survival of binary neutron star systems may in part be due to their lower kick speeds. We also find a correlation between the kick and both the mass dipole of the ejecta and the explosion energy. Furthermore, one channel of black hole birth leaves masses of ∼10M_⊙, is not accompanied by a neutrino-driven explosion, and experiences small kicks. A second channel is through a vigorous explosion that leaves behind a black hole with a mass of ∼3.0M_⊙kicked to high speeds. We find that the induced spins of nascent neutron stars range from seconds to ∼10 ms, but do not yet see a significant spin/kick correlation for pulsars. We suggest that if an initial spin biases the explosion direction, a spin/kick correlation would be a common byproduct of the neutrino mechanism of core-collapse supernovae. Finally, the induced spin in explosive black hole formation is likely large and in the collapsar range. This new 3D model suite provides a greatly expanded perspective and appears to explain some observed pulsar properties by default.

We present in this paper a public data release of an unprecedentedly large set of core-collapse supernova (CCSN) neutrino emission models, comprising 100 detailed 2D axisymmetric radiation-hydrodynamic simulations evolved out to as late as ∼5 s post-bounce and spanning an extensive range of massive-star progenitors. The motivation for this paper is to provide a physically and numerically uniform benchmark data set to the broader neutrino detection community to help it characterize and optimize subsurface facilities for what is likely to be a once-in-a-lifetime galactic supernova burst event. With this release, we hope to (1) help the international experiment and modelling communities more efficiently optimize the retrieval of physical information about the next galactic CCSN, (2) facilitate the better understanding of core-collapse theory and modelling among interested experimentalists, and (3) help further integrate the broader supernova neutrino community.

We have simulated the collapse and evolution of the core of a solar-metallicity 40M_⊙star and find that it explodes vigorously by the neutrino mechanism, despite its very high “compactness.” Within ∼1.5 s of explosion, a black hole forms. The explosion is very asymmetrical and has a total explosion energy of ∼1.6 × 10⁵¹erg. At black hole formation, its baryon mass is ∼2.434M_⊙and gravitational mass is 2.286M_⊙. Seven seconds after black hole formation, an additional ∼0.2M_⊙is accreted, leaving a black hole baryon mass of ∼2.63M_⊙. A disk forms around the proto−neutron star, from which a pair of neutrino-driven jets emanates. These jets accelerate some of the matter up to speeds of ∼45,000 km s⁻¹and contain matter with entropies of ∼50. The large spatial asymmetry in the explosion results in a residual black hole recoil speed of ∼1000 km s⁻¹. This novel black hole formation channel now joins the other black hole formation channel between ∼12 and ∼15M_⊙discovered previously and implies that the black hole/neutron star birth ratio for solar-metallicity stars could be ∼20%. However, one channel leaves black holes in perhaps the ∼5–15M_⊙range with low kick speeds, while the other leaves black holes in perhaps the ∼2.5–3.0M_⊙mass range with high kick speeds. However, even ∼8.8 s after core bounce the newly formed black hole is still accreting at a rate of ∼2 × 10⁻²M_⊙s⁻¹, and whether the black hole eventually achieves a significantly larger mass over time is yet to be determined.

In this paper, we analyze the neutrino-driven winds that emerge in 12 unprecedentedly long-duration 3D core-collapse supernova simulations done using the code Fornax. The 12 models cover progenitors with zero-age main-sequence mass between 9 and 60 solar masses. In all our models, we see transonic outflows that are at least 2 times as fast as the surrounding ejecta and that originate generically from a proto−neutron star surface atmosphere that is turbulent and rotating. We find that winds are common features of 3D simulations, even if there is anisotropic early infall. We find that the basic dynamical properties of 3D winds behave qualitatively similarly to those inferred in the past using simpler 1D models, but that the shape of the emergent wind can be deformed, very aspherical, and channeled by its environment. The thermal properties of winds for less massive progenitors very approximately recapitulate the 1D stationary solutions, while for more massive progenitors they deviate significantly owing to aspherical accretion. TheY_etemporal evolution in winds is stochastic, and there can be some neutron-rich phases. Though no strongr-process is seen in any model, a weakr-process can be produced, and isotopes up to⁹⁰Zr are synthesized in some models. Finally, we find that there is at most a few percent of a solar mass in the integrated wind component, while the energy carried by the wind itself can be as much as 10%–20% of the total explosion energy.

We present updated atmospheric tables suitable for calculating the post-formation evolution and cooling of Jupiter and Jupiter-like exoplanets. These tables are generated using a 1D radiative transfer modeling code that incorporates the latest opacities and realistic prescriptions for stellar irradiation and ammonia clouds. To ensure the accuracy of our model parameters, we calibrate them against the measured temperature structure and geometric albedo spectrum of Jupiter, its effective temperature, and its inferred internal temperature. As a test case, we calculate the cooling history of Jupiter using an adiabatic and homogeneous interior and compare with extant models now used to evolve Jupiter and the giant planets. We find that our model reasonably matches Jupiter after evolving a hot-start initial condition to the present age of the solar system, with a discrepancy in brightness temperature/radius within 2%. Our algorithm allows us to customize for different cloud, irradiation, and metallicity parameters. This class of boundary conditions can be used to study the evolution of solar system giant planets and exoplanets with more complicated interior structures and nonadiabatic, inhomogeneous internal profiles.

The two-moment method is widely used to approximate the full neutrino transport equation in core-collapse supernova (CCSN) simulations, and different closures lead to subtle differences in the simulation results. In this paper, we compare the effects of closure choices on various physical quantities in 1D and 2D time-dependent CCSN simulations with our multigroup radiation hydrodynamics code Fornax. We find that choices of the third-order closure relations influence the time-dependent simulations only slightly. Choices of the second-order closure relation have larger consequences than choices of the third-order closure, but these are still small compared to the remaining variations due to ambiguities in some physical inputs such as the nuclear equation of state. We also find that deviations in Eddington factors are not monotonically related to deviations in physical quantities, which means that simply comparing the Eddington factors does not inform one concerning which closure is better.